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Abstract. 

The NASA Orbital Debris Observatory (NODO) astronomical survey uses 
a transit 3-m liquid mirror telescope to observe a strip of sky in 20 
narrow-band filters. In this article, we analyze a subset of data from the 
1996 observing season. The catalog consists of 18,000 objects with 
10<V<19 observed in 10 narrow-band filters ranging from 500 nm to 950 
nm. We first demonstrate the reliability of the data by fitting the 
Bahcall-Soneira model of the Galaxy to the NODO magnitude counts 
and color counts at various galactic latitudes. We then perform a 
hierarchical clustering analysis on the sample to extract 206 objects, out 
of a total of 18,000, showing peculiar spectral energy distributions. It is a 
measure of the reliability of the instrument that we extract so few 
peculiar objects. Although the data and results, per se, may not seem 
otherwise particularly remarkable, this work constitutes a milestone in 
optical astronomy since this is the first article that demonstrates 
astronomical research with a radically new type of mirror. 

Subject headings: 

instrumentation: miscellaneous, surveys, stars: peculiar, Galaxy: general, 
galaxies: peculiar, quasars: general. 
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1. INTRODUCTION 



It has been known for centuries that the surface of a liquid 
spinning in a gravitational field has a parabolic shape. In principle, a 
spinning container filled with a reflecting liquid such as mercury could 
then be used as the primary mirror of a telescope. In practice, the concept 
never was taken seriously by astronomers for a variety of reasons, 
foremost of which are the facts that past experiments met with 
difficulties and, especially, that the mirror cannot be tilted and hence 
cannot track like a glass mirror. In recent times, the concept has been 
resurrected by Borra (1982) who argued that modern technology renders 
liquid mirrors useful to astronomy. Since that article, several 
developments have shown that the concept is viable. Optical tests 
showed that a 1.5-m diameter liquid mirror (LMT) was diffraction- 
limited (Borra et al. 1992): This article gives a wealth of information on 
the basic LM technology. This was followed by tests of a 2.5-m diameter 
liquid mirror (Girard & Borra (1997)) that also showed diffraction limited 
performance. Independent optical tests of a liquid mirror by Ninane & 
Jamar (1996) have confirmed diffraction-limited performance. Deep-sky 
imagery, obtained with the UBC-Laval 2.7-m LMT, has been 
demonstrated by Hickson et al. (1994). 

Liquid mirrors are interesting in other areas of science 
besides Astronomy. For example, the University of Western Ontario 
has built a Lidar facility that houses a 2.65-m diameter liquid mirror as 
receiver (Sica et al. 1995). This telescope has now been operated trouble 
free for 5 years in the harsh Canadian climate, demonstrating the 
robustness of LMTs. The telescope is however not an imager. A Lidar 
facility has also been built and operated in Alaska by the University of 
California at Los Angeles. The most successful imager is the NODO 3-m 
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diameter LMT (http://www.sunspot.noao.edu/Nodo/nodo.html) built 
by NASA (Potter & Mulrooney 1997, Mulrooney 1998) to monitor space 
debris. The observatory has operated for a few years. NODO has also 
been used to obtain a large number of nights of astronomical data 
(Hickson & Mulrooney 1997). Liquid mirrors have also interesting 
engineering applications (Thibault, Szapiel & Borra, 1998). 

Although now somewhat outdated in this rapidly evolving 
field, a review paper (Borra 1995) elaborates on many aspects of Liquid 
Mirrors and Liquid Mirror Telescopes (including the all-important issue 
of the field of regard). Useful reference papers on LM technology and 
related issues have been written by Hickson and collaborators, ( e.g. 
Hickson et al. (1993), Gibson & Hickson (1992a, 1992b)). 

While optical shop tests are needed to evaluate the intrinsic 
quality of the mirror and pretty pictures are reassuring, members of the 
astronomical community have repeatedly emphasized through personal 
communications that "the proof of the pudding is in the eating"; hence 
at some stage it is necessary to demonstrate scientific research. A 
previous article (Content et al. 1989) published astronomical research 
with a LMT but it used a small make-shift telescope with rudimentary 
instrumentation. The present paper demonstrates astronomical research 
with data obtained with a professional-quality 3-m LMT and is another 
important milestone in the demonstration of LMTs as Astronomical 
instruments. 



2. DATA DESCRIPTION 

The primary data set consists of 34 nights of narrow-band 
observations taken during the spring of 1996. Each night uses a single 
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narrow-band filter selected among 11 bands (figure 1). The observations 
were taken with the NASA Orbital Debris Observatory (NODO) and are a 
subset of those described by Hickson & Mulrooney (1997). However, 
although we used the same primary images, we did our own 
independent data reduction with different algorithms and software, as 
described in the next section. The NASA Orbital Debris Observatory uses 
a 3-m diameter transit LMT. It is equipped with a 2048x2048 thick CCD 
camera continuously read-out at sidereal speed in the Time Delayed 
Integration mode (TDI). Tracking is thus done electronically by stepping, 
on the CCD, the pixels in the East-West direction at the sidereal rate. The 
integration time per night per object (97 seconds) is given by the time it 
takes an object to cross the 20-arcmin field of view. Longer integration 
times could be achieved by coadding data from different nights; 
however, in this work, we did not do any coadding. A typical night of 
observation gives a continuous band of observation which is divided, for 
convenience, in 280 20X20 arcminutes frames at constant declination, 
and increasing right ascension for a total of about 30 square degrees. The 
declination is given by the latitude of the observatory (32.96 deg). Each 
frame takes 8 Mb for a total of 2.2 Gb per night of observation. 

The primary data has been described by Hickson & Mulrooney 
(1997) and evaluated by them so that there is no need to repeat their 
work here. Figure 2 gives an overview of the 27 nights of data that we 
examined and illustrates the robustness of the system. Every night of 
observation is labeled with the central wavelength of the bandpass used. 
The abscissa gives the time of observation converted to right ascension 
(top) as well as galactic latitude and galactic longitude (bottom). The 
timeline of a given night of observation is colored white whenever the 
images were deemed to be usable. The line is colored gray whenever the 
data was found to be unusable for whatever reason, including clouds. At 
the right of the timeline a brief comment is made. We can readily see 
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that the majority of the timelines are white. The gray-shaded portions of 
the timelines are usually due to either clouds, to twilight or to cooling 
problems with the CCD detector which was too warm, hence had 
excessive dark counts. Although most of the nights were reasonably 
clear, some of them were definitely not photometric. In order to limit 
possible errors in the photometry from the beginning to the end of a 
night, we chose to extract the information from the best 11 nights, one 
per filter. Figure 1 gives the transmission curve of the filters , the 
quantum efficiency of the CCD (Hickson & Mulrooney 1997), and the 
night sky emission at Kitt Peak and Observatoire de Haute-Provence 
(Massey et al. 1990). 

3. DATA REDUCTION 

The entire data set was roughly 60 Terabytes in size. In order to 
extract useful information from this huge amount of data we wrote a set 
of C routines optimized for the format of the images. Each night was cut 
into about 280-300 2kx2k frames with an overlap of one arcmin between 
them. The set of automated C routines reduced the data according to the 
following steps: 

3.1 Image analysis 

0. Preprocessing was already performed on our images by Hickson & 
Mulrooney (1997) who removed the dark counts and subtracted the 
bias. We did not flat field our frames since driftscanning gives a very 
flat response. In the East- West direction, variations are near zero and 
we would not expect variations of more then 1% in addition to the 
gradient in the North-South direction. 
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1. Remove bad pixels and bad columns and average the gap using the 
nearest pixels and columns. This step is quite straightforward once 
one knows the coordinates of the faulty pixels. The CCD matrix was 
good except for one two-pixel column. 

2. Convolve the image with a 3x3 pix gaussian kernel (FWHM = 2.35 pix 
= 1.4 arcsec). This reduces sky noise when present and increases the 
number of bits in the sky histogram. Sky flux was not always available 
because the CCD was not sensitive enough in the blue (fig. 1) and/ or 
the integration time was too small. As a result, the noise was not sky- 
dominated for the blue filters (X<6000A) but rather electronics- 
dominated. This did not have dramatic effects in the detection process, 
but it prevented us from measuring the actual magnitude limit of the 
sky and total efficiency of the system. Sky values and standard 
deviations are then measured and computed in a mosaic of 32x32 
elements, each of them 64x64-pix large. The size of the mosaic 
elements is optimized to get the finest possible sky value within a 
region where the large-scale sky variation is lower then 1 a. 

3. Detect all pixels above 2 a of sky value (from the 32x32 mosaic) and 
identify groups bigger then 9 contiguous pixels and tag them as objects. 
This 2-g threshold, is our primary criterion to distinguish real objects 
from false detections. The 9 pixels can form any shape as long as they 
are contiguous. Hence, the total signal-to-noise ratio of the faintest 
detection is 6. This threshold is actually low since it allows 50% to 80% 
of all detections to be false. False detections are later eliminated by 
merging the nights. 

4. Measure centroid, axial ratio, angle, Kron radius and growth curve 
from 2 to 10 pixels around every objects. The centroid is computed 
using the center of gravity of the complete distribution of all 
contiguous pixels. The axial ratio is calculated from all contiguous 
pixels without a priori subtraction. The position angle is the angle 



7 



between the semi-major axis and the north-south direction. The axial 
ratio and the position angle are computed using the second moment 
of the spatial distribution (Hickson & Mulrooney 1997). The Kron 
radius is the isophotal radius defined in Kron (1980). The growth 
curve was computed over a circular disk between r = 2 pix and r = 10 
pix by steps of 1 pix. Most of the objects are stars, thus point-like, and 
the axial ratio and position angles are usually sensitive to the wings of 
the PSFs near the edge of the field where the TDI sidereal velocity is 
either too small (north) or too large (south). Where Hickson & 
Mulrooney (1997) use a corrected isophotal flux, we use a simple 
aperture photometry with a 3 arcsec diameter circular aperture. 

3.2 Photometric and astrometric calibration 

The photometric characteristics of the data, albeit reduced with 
different software and algorithms, have been described by Hickson & 
Mulrooney (1997) and evaluated by them; therefore we will not repeat it 
here. We shall, however, describe briefly our own data reductions. 

To make a zero-point calibration of the nights we searched 
astronomical databases to find calibrated standard stars in our strip of 
sky. We only found two blue stars, one at each end of the strip, with 
featureless spectra calibrated in AB magnitudes (Oke 1990). In practice, 
only the star at the beginning of the night could be used since the one at 
the end of the night was too bright and saturated the CCD in some filters. 
Tests show that the zero-point calibration does not vary by more the 0.3 
mag from the beginning to the end of the nights that we chose. 
Fortunately, the search for peculiar objects, which is the primary 
scientific objective of this article does not need an accurate zero-point 
calibration. 
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After the primary reduction and photometric measurements were 
done, the next important step was to carry out an accurate astrometry 
over the whole night. The coordinates of the objects are the principal 
criterion used for identification and merging of objects from one night to 
another. The astrometry is calibrated in two steps, first we use analytic 
transformations from instrumental coordinates to J2000 equatorial 
coordinates. This allows to correct for precession, nutation, and 
aberration. Because LMTs are zenith instruments, they are particularly 
sensitive to astronomical aberration (up to 20 arcsec from the beginning 
to the middle of the night). For the second step we wrote a fine-tuning 
program, using standards from the United States Naval Observatory 
astrometric survey (Monet 1997). We corrected the coordinates frame by 
frame comparing all matched stars (usually over 50 /frame, up to 500). 
The final astrometry is good within 0.5 arcsec. We merge the nights into 
one catalog containing the coordinates and the magnitudes of all objects. 
We only included in the catalog the objects detected in every filter. 
Figure 3 shows the magnitude counts in all of the filters of the final 
catalog. It is difficult to estimate the completeness of the catalog because 
the quantum efficiency of the detector varies considerably over the 4,000 
A spectrum. In particular, the CCD is not very sensitive in the blue; 
hence the catalog is biased toward faint blue objects. 

4. REPRODUCING THE KNOWN: COMPARISONS TO IMAGES FROM 
THE PALOMAR OBSERVATORY SKY SURVEY AND TO MODELS OF 
THE GALAXY 

Before a new instrument is used to carry out original research it 
must demonstrate that it can reproduce known results. We therefore 
first make a qualitative comparison to images from the Digitized Sky 
Survey II. We then compare our star counts with standard models of the 
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Galaxy. Hickson & Mulrooney (1997) also give an evaluation of the data 
which complements ours. 

4.1 Comparing the LMT images to images from the Palomar Observatory 
Sky Survey. 

Eye inspections of a very large number of frames in all filters and 
at all times of the nights indicate that the images are well behaved. We 
also visually compared in greater details selected frames to the 
corresponding images of the DSS II. For example, figure 4 shows such a 
comparison between an unfiltered (equivalent to a broad-band red filter) 
deep image of the NODO survey (top) and the DSS II image (bottom). 
One may appreciate the similarity between the CCD image and the 
digitized photographic plate. NODO narrow-band images are generally 
deeper in red than DSS II, and shallower in blue. Hickson & Mulrooney 
(1997) gives limiting magnitudes and associated magnitude errors for all 
filters. 

4.2 Comparing star counts and colors to models of the galaxy 

We next compare star counts, extracted from our database to the 
well-established Bahcall-Soneira models of the Galaxy. Bahcall (1986) 
proposed a simple analytical model of the Galaxy containing a spheroid 
and a disk component and wrote a code computing the magnitude 
counts and number counts by simply summing the contributing stellar 
populations in the line-of-sight of the observer. The Bahcall-Soneira 
models reproduce past observations remarkably well (Bahcall et al. 1985, 
Lasker et al. 1987), and they can be considered as reliable standards for the 
distribution of the stars in the Galaxy in a range of magnitude B < 20. 
Recent studies by Santiago et al. (1996) using HST data, Bath et al. (1996) 
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using the APS catalog, and Malkov & Smirnov (1994) using the Guide 
Star Catalog show similar agreements in star counts, but some 
discrepancies in color counts attributed to a thick disk component. We 
used the Galaxy parameters given in Table 1 (Bahcall 1986). 

The code was found at http:/ /www.sns.ias.edu/jnb/ galaxy/html, 
courtesy of J. Bahcall. We did not change the normalization of its 
computed results. We use directly the output of the program to produce 
count and color curves at three different galactic latitudes for different 
magnitude intervals. We transform the standard B and V magnitudes to 
AB magnitudes in our narrow-band filters at 500 and 700 nm. The 
transformation equations are calculated using the standard spectra of the 
blue star Hz 21 from Oke (1990) for the zero-point calibration, and the 
Gunn & Stryker catalog of stellar spectra (1983) for the color calibration. 
Least square fits yield the transformation equations: 



M50 = V + 0.385 (B-V) - 0.033, 
and for -0.4<B-V<0.9: 
(M50-M75) = 1.36 (B-V) - 0.40 
for B-V>-0.9: 

(M50-M75) = 3.00 (B-V) - 2.1 



(1) 
(2) 
(3) 



where M50 and M75 are the AB magnitudes in our catalog, and B and V 
are the standard bands used by Bahcall and Soneira. Figures 5 and 6 show 
the data and models for three different galactic latitudes and longitudes 
where extinction is less than 0.01 (Burstein & Heiles 1982). In the models, 
the counts are sensitive to the giant branch population of the spheroid 
component of the Galaxy. The metal poor cluster M15 gives the best fits 
towards the north galactic pole, and the rich cluster 47 Tucanae gives the 
best fits for lower galactic latitudes. This is what one expects if spheroid 
stars are more metal-poor than disk stars. Our data is thus consistent 
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with past observations. The models fit well the counts for four different 
magnitude intervals, although photometric errors make the comparison 
less reliable for V>19. There is a small shift in the color peak near the 
plane of the galaxy (corresponding to the end of the nights). Our data are 
plotted as they come straight out of the database. There is no 
renormalization in neither counts nor colors. Hence, because we only 
use a single photometric standard at the beginning of the nights, the shift 
is probably due to extinction variations from the beginning to the end of 
the nights. It is a testimony to the photometric qualities of the nights that 
the comparisons between models and data are as good as shown in 
figures 5 & 6. 

5. A SEARCH FOR PECULIAR OBJECTS WITH A HIERARCHICAL 
CLUSTERING ANALYSIS 

5.1 Simulations 

We created a database of artificial spectral energy distributions of 
different kinds of astronomical objects, as they would look like seen 
through our filters, by convolving higher resolution spectra, obtained 
from the literature, with the 10 filters shown in figure 1. For stars, we 
used the Gunn & Stryker catalog of stellar spectra . For the quasar we 
used a composite "typical" spectrum. Unfortunately quasar spectra 
display a broad variety of spectra and a "typical" spectrum is a poor 
approximation. Figure 7 shows our library spectra for E , Sc and Irr 
galaxies at the different redshifts as well as our "typical" QSO at different 
redshifts. The galactic spectra have not been corrected for evolution. We 
assume that whatever might happen in the stellar population of the 
giant galaxies, it will only change the continuum by small amounts for 
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z<0.5. This is probably a crude approximation for dwarf irregulars and 
galaxies that experience starbursts, but this assumption is not critical, 
since the purpose of the simulation is not to test evolutionary models 
but to get an idea of what galaxy spectra, seen through our filters, would 
look like at different redshifts. Figure 7 shows that Irregular galaxies and 
Sc galaxies as well as QSOs have striking features especially at high 
redshift. Since the filters do not cover the whole range between 500nm to 
950nm, emission lines are not always detected. For instance, while the 
QSO shows strong lines at z=3.9, because both Ly a (1216A) and CIV 
(1549A) fall in two of our filters, the same lines are no longer visible at 

o 

z=4.1. Similarly, the usually conspicuous break at 4000 A in giant 
ellipticals (Cabanac et al. 1995) is almost invisible at z=0.3, in figure 7. The 
dwarf irregular galaxies show an interesting continuum inversion at 
high redshift (z>0.45). This simulation gives another tool to extract 
information from the peculiar object spectra given by the HCA (next 
section). 

5.2 Hierarchical Clustering Analysis 

In the rest of the paper we will refer to flux distributions of the 10 
filters as spectra, although we do not uniformly cover the entire range 



between 4000 A and 9000 A (see figure 1) and thus do not have the 
complete low-resolution spectra of the objects. Notwithstanding these 
limitations, our catalog carries more information than an equivalent 
BVRI photometry catalog. Here, we present a search for peculiar objects 
on a subset of 18,000 bright objects using a hierarchical clustering analysis 
(hereafter HCA). The method is based on the technique used by 
Beauchemin et al. (1993) (also Beauchemin et al. 1991) using HCA code 
from Murtagh & Heck (1987). A comprehensive description of the 
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technique can be found in the previous cited papers and references 
therein. HCA is a multi-dimensional clustering technique that uses a 
minimum variance criterion to segregate objects in groups. To perform 
HCA, we sliced the catalog into sub-samples of 2,000 objects each. Table 2 
gives the beginning and the ending right ascensions, and total area for 
each sub-sample. The total area of each sub-sample decreases from the 
north galactic pole to the galactic plane as a result of an increasing 
density of objects. We sliced our catalog for two reasons: First, HCA is 
time-consuming and the CPU cost increases as the square of the number 
of objects. Second, each sub-sample covers a shorter time length then the 
full night, slicing the data set reduces the impact of photometric 
variations. Our analysis indicates that the impact of varying atmospheric 
extinction is negligible within our time subsamples. We excluded the 
flux of the last filter (#95) from the objects, because the combined low 
sensitivity of the CCD and high sky brightness (figure 1) induced large 
photometric errors which falsely produced many peculiar objects. We 
did not select any stopping rule to choose the optimal number of clusters 
(Milligan et al. 1985) since our purpose was not to create a rigorous 
classification of objects but rather to group them according to a relative 
spectral morphology within a given number of clusters. The HCA 
output gives a hierarchical tree (or dendogram) with a maximum of 40 
branches (or clusters) assembling similar objects. We may cut the 
dendogram at any of the 40 levels, corresponding to the number of 
clusters we want. To choose the number of clusters, we tested the 
dendograms with 5, 10, 20 and 40 clusters. We needed to isolate as many 
peculiar objects as possible from the rest of the subset. We concluded that 
the best number of clusters was 40. A smaller number of clusters did not 
leave enough free clusters for outliers. Figure 8 gives the result of HCA 
on the group 0. Median spectra of selected groups are plotted in figure 9. 
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For the HCA, the objects are normalized at 700 nm. We plot the median 
spectra with their true AB magnitudes. 

All objects in clusters containing less than 10 objects were tagged as 
potentially interesting candidates. In this first iteration, we extracted 445 
objects from the 18,000 objects of the catalog. We applied the HCA a 
second time on the sub-sample, using 10 clusters. This effectively 
classified the sub-sample by colors. Figure 10 shows the result of the HCA 
and median spectra of each cluster are shown figure 11. Figure 12 
compares 5 observed stellar spectra to the simulations. The observed 
spectra reproduce well the simulated spectra, another testimony as to the 
reliability of our data. 

In order to associate an approximate spectral type to each object, we 
fitted the 445 objects with the best fitted spectrum from a database of 
spectral energy distributions obtained from our spectral library (section 
5.1). Even if our spectra were not rigorously corrected, this step 
permitted to isolate all spectra that looked stellar. We eliminated from 
our first sample of 445 objects all objects showing unambiguous M and K 
star features. Because cool stars are relatively rare at these magnitudes, 
the HCA classified a number of these normal stars as peculiar. 

Because different wavelengths are observed on different nights, 
variable objects have unreliable spectra. Our spectra are also light curves! 
Note however that consecutive filters were not taken on consecutive 
nights so that time does not increase (nor decrease) from blue to red. 
Light variations will show up as spectral variations so that we should 
expect that a fraction of our peculiar objects are actually variable objects. 
This is almost certainly the case for stellar looking continua that show a 
deep absorption (see figure 13, #126, #38). Objects that have strong 
enough light variations show up as peculiar spectra and our catalog is 
probably heavily contaminated by these objects. We remove from our 
sample some of the most obvious cases. 
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Finally, we keep the remaining 206 objects and looked at their 
spectra and images for hints of peculiar features, emission peaks, halos, 
or companions. The list of these peculiar objects with their coordinates, 
AB magnitudes in filter 550 (essentially the V band) comments such as, 
suspected type, best-fitted stellar library spectra, striking spectral features 
is given in Table 3. A postscript figure containing all peculiar spectra is 
available at 

http:/ /wood.phy.ulaval.ca/lmt/ abstracts/peculiar.spectra.ps.z (gzipped 
230kb). Figure 13 illustrates what we mean by absorption lines (a) or (sa), 
emission lines (e) or (se), breaks (b) or (sb), variable (v), and inverted 
continuum (i) in the comments of table 3 ("s" in sa, se and sb stands for 
"strong"). Eyeball classification of spectral features is highly subjective 
since it depends on where the eye sets the continuum. Independent 
classification by two of us shows disagreements as to whether a feature is 
an absorption, an emission or a break. Consequently, the classification of 
a feature in Table 3 should be taken cautiously: the comments really only 
indicates that there is some spectral feature near that wavelength. 
Follow-up spectroscopy is needed. The following section is devoted to 
particularly interesting objects that deserve a spectroscopic follow-up. 

5.3 Peculiar Objects 

In this subsection, we comment on some of our most interesting 
objects. 

Optical counterparts of radio source : Objects #50, 122 

We checked whether some of our 206 peculiar objects coincide 
with known peculiar objects by cross-referencing them with the 
NASA/IP AC Extragalactic Database (NED). We find that two of them are 
optical counterparts of known radio sources. Figure 14 shows the red 
image of each object with superimposed contours of the Radio Source. 
The image comes from the Digital Sky Survey and the Radio contours 
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come from the VLA FIRST Survey. Both optical images and radio 
contours are extracted from the SkyView image database. The radio 
sources, 87GB123819.7+332002 (obj#50; already known as an optical 
counterpart in NED) and 87GB1 62630.1 +325701 (obj#122; discovered by 
us) have steep radio spectral indices, implying non-thermal emission 
processes in the transparent radio lobes. Both optical objects show 
continuum radio emission at 1.4 Ghz. Extended emissions in the 
immediate vicinity of the objects are barely detectable on the DSS. The 
south lobe of the obj#50 seems to have an optical counterpart, and the 
obj#122 (galaxy) is surrounded by what looks like faint galaxies. Obj#122 
might be of particular interest since the lobes have peculiar 
morphologies and dynamical studies could infer properties of the galaxy 
cluster. 

Gravitational lens candidates , close pairs and groups: Objects #48-49, 
176-177 

We find two gravitational lens candidates(#48-49, 176-177), as well 
as several close pairs and groups are identified (#53-54, 40-119, 57-61-146, 
153-157, 169-170, 129-134). The selection criteria are that the objects 
should have similar spectra and that their centroids should be separated 
by less than 10 arcsec for gravitational lens candidates, and 10 arcmin for 
close groups. Figure 15 plots the spectra of the two gravitational lens 
candidates. The close groups are probably stars with similar spectra, but 
they could also be remote members of galaxy clusters. At a redshift of 1, 
an angular distance of 10 arcmin is equivalent to a distance of 5 Mpc 
(q =l/2, H^SOkm.s^Mpc 1 ). 
Blue objects : Objects #1-47 

14 blue objects (#1-14) were detected with the colors of O or B stars, 
and 33 objects (#15-47, 185) with colors of A stars. They could be hot 
subdwarfs, white dwarfs or QSOs. 
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Object #203: Galaxy cluster 

Finally, an interesting object was found while doing eyeball 
checking of a randomly selected sample of Table 3. Obj#203 is a red 
object (probably an elliptical galaxy) surrounded by a great number of 
faint extended objects. The DSS image reveals a large number of galaxies 
inside a radius of 10 arcmin around obj#203. 
Inverted spectrum objects 

These objects (e.g., figure 13, #150) are potentially quite interesting; 
however we find that the majority of them are very faint, so that the 
inversion may actually come from faulty sky subtractions. However, at 
least two (#14, #134) are bright enough that this explanation does not 
hold. We cannot exclude that they are variable objects; however the fact 
that increasing time does not correlate with increasing (nor decreasing) 
wavelengths, makes this somewhat unlikely. 

6. CONCLUSION 

The main scientific goal of this paper consists in a search for 
objects that have peculiar energy distributions. We extract spectra for 
18,000 objects having 10<V<19 in a 15-arcminute strip of sky extending 
from 12h30 to 18h00 for a total of -18 square degrees. We then perform a 
hierarchical clustering analysis on the sample in order to extract objects 
showing peculiar spectral energy distributions. A blind HCA analysis 
finds 445 objects having peculiar energy distributions. Inspection of the 
energy distributions shows that some of these objects are actually cool 
stars so that they were removed from the list. Because different filters are 
used on different nights, variable objects show up as peculiar objects. We 
therefore removed from the list objects that are obvious variable stars as 
well as some objects that have false features caused by problems with the 
detector (e.g. cosmic rays). The final list contains 206 objects. Because we 
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do not have redundant observations, there is a very high likelihood that 
most of these objects are actually variable objects, so that follow-up 
spectroscopy is needed to confirm their peculiarities. On the other hand, 
to the best of our knowledge, this survey is unique because of the 
combination of spectral resolution, relatively faint limits and area 
surveyed. It therefore has good potential to find rare objects. 

Although the data and results, per se, may not seem particularly 
remarkable, this paper constitutes a milestone in optical astronomy since 
it is the first paper to demonstrate astronomical research with such a 
radically new type of optics. This is a first generation liquid mirror 
telescope and its instrumentation is actually not optimized for 
astronomical observations since the detector is a thick CCD with low 
quantum efficiency and the NODO 3-m telescope was designed to 
observe fast moving space debris so that the electronics generates a 
strong read-out noise. Still, the telescope demonstrated its robustness by 
consistently giving good images over a large number of nights. That the 
data is reliable is demonstrated by the fact that star counts and color 
histograms fit the Bahcall & Soneira model of the Galaxy. It is also 
demonstrated by the fact that the impartial HCA search only found 445 
objects out of 18,000 in its first blind pass. 

Considering the technology at the time of this writing, the low 
capital and operational costs of the first generation liquid mirror 
telescopes make them ideal specialized instruments for narrowly focused 
surveys where the limited field of regard is not a serious limitation. 
However, future advances in corrector designs could render the next 
generation LMTs far more versatile by greatly extending their fields of 
regard (Borra, 1995). 
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Figure captions 
Figure 1 

It gives the typical quantum efficiency (QE) for a thick 2kx2k, the sky 
brightness from Kitt Peak, and the transmission curves (TC) of the filters. 

Figure 2 

Overview of a subset of the 1996 season of observation of the NASA 
Orbital Debris Observatory. The gray shades overlap non-optimal 
conditions of observation. 

Figure 3 

Plots the magnitude counts of the final catalog for objects detected in all 
colors. The completeness is poor beyond mag~18.5 in red and 17.5 in 
blue. This comes from the CCD which is not optimized for astronomical 
purposes. There are combined effects of saturation in red (mag<13) and 
poor efficiency in blue that strongly limit the linear range of detection. 

Figure 4 

Comparison between an image taken with the NASA Orbital Debris 
Observatory (top) and its POSS II counterpart (bottom). The field is 5 x 7 
arcmin. RA = 12h08m, DEC = 33.0° Q2000.0). 

Figure 5 

Plots magnitude counts at three galactic latitudes: (1=86°, b=83°) (1=53°, 
b=47°) (1=56°, b=34°). The steps are the observed counts per square 
degree, and the dash-lines are the Bahcall-Soneira models. 

Figure 6 

Plots the observed color counts (in steps) at same latitudes than figure 5, 
for four intervals of magnitudes (V<14, 14<V<17, 17<V<18.5, 
18.5<V<19.5) and their associated Bahcall-Soneira models of the Galaxy 
(in dash-lines). The different theoretical curves use different stellar 
populations for the spheroid component of the Galaxy. 

Figure 7 

Gives the result of HCA on group (first 2000 objects) of the catalog. The 
histogram shows the number of objects in each cluster. 

Figure 8 

Plots of the median spectra for each cluster of the Fig. 7. 
Figure 9 
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Shows the result of HCA for the 445 objects extracted from the 18,000 
objects. 

Figure 10 

The plots show the median spectra for each of the 10 clusters of figure 8. 
Figure 11 

Show simulations of stellar spectra (B, F, G, K, M) filtered in our 10 
bands, different types of galaxies (E, Sc, Irr) at various redshifts and a 
composite quasi stellar object spectra at various redshifts. It illustrates, 
the diversity of the possible objects (no evolution). 

Figure 12 

Gives 5 spectra and their best-fitted stellar spectra, from 08 to M8. All 
objects which showed obvious similarities with stars (e.g. M8) were 
removed from the peculiar object catalog. 

Figure 13 

Illustrates the comments given in Table 3. Absorption (a or sa), emission 
(e or se), break (b or sb), inverted continuum (i) and variable (v). the 
numbers are the same as in Table 3. 

Figure 14 

Images of two of the peculiar objects which are optical counterpart (line- 
of-sight) of two extragalactic radio sources. The gray scale images are 
DSSII images and the contours are FIRST 1.4Ghz flux. 

Figure 15 

Spectra of two gravitational lens candidates. 
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Table titles 
Table 1: 

Galaxy parameters for the model of the galaxy. 
Table 2: 

HCA sub-samples. 
Table 3: 

Peculiar objects 
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